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The magnetization, M(H < 30 T, 0.7 K < T < 300 K), of (C5Hi2N)2CuBr4 has been used to 
identify this system as an S = 1/2 Heisenberg two-leg ladder in the strong-coupling limit, J± — 13.3 
K and Jy = 3.8 K, with Hci = 6.6 T and Hc2 = 14.6 T. An inflection point in M{H,T = 0.7 K) at 
half-saturation, Ms/2, is described by an effective XXZ chain. The data exhibit universal scaling 
behavior in the vicinity of Hci and Hc2, indicating the system is near a quantum critical point. 
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Magnetic spin ladders are a class of low dimensional 
materials with structural and physical properties between 
those of ID chains and 2D planes. In a spin ladder, the 
vertices possess unpaired spins that interact along the 
legs via J|| and along the rungs via J±, but are isolated 
from equivalent sites on adjacent ladders, i.e. interlad- 
der J' <^ J\\,J±- Recently, a considerable amount of 
attention has been given to the theoretical and experi- 
mental investigation of spin ladder systems as a result 
of the observation that the microscopic mechanisms in 
these systems may be related to the ones governing high 
temperature superconductivity [|l]j2) . The phase diagram 
of the antiferromagnetic spin ladder in the presence of a 
magnetic field is particularly interesting. At T = with 
no external applied field, the ground state is a gapped, 
disordered quantum spin liquid. At a field Hd, there is 
a transition to a gapless Luttingcr liquid phase, with a 
further transition at Hc2 to a fully polarized state. Both 
Hci and Hc2 are quantum critical points [||. Near these 
points, the magnetization has been predicted to obey a 
universal scaling function but until now, this behavior 
has not been observed experimentally. 

A number of solid-state materials have been pro- 
posed as examples of spin ladder systems, and an ex- 
tensive set of experiments have been performed on 
Cu2(C5Hi2N2)2Cl4, referred to as Cu(Hp)Cl §. The 
initial work identified this material as a two- leg 5=1/2 
spin ladder |4|-[l2|. Although quantum critical behav- 
ior has been preliminarily identified in this system near 
Hci , this assertion is based on the use of scaling parame- 
ters identified from the experimental data rather than the 
ones predicted theoretically ||ll|,f2|. Furthermore, more 



recent work has debated the appropriate classification of 
the low temperature properties |13-^]. Clearly, addi- 
tional physical systems are necessary to experimentally 
test the predictions of the various theoretical treatments 
of two-leg S = 1/2 spin ladders. 

Herein, we report evidence that identifies bis(piperi- 
dinium)tctrabromocuprate(II), (C5Hi2N)2CuBr4 
hereafter referred to as BPCB, as a two-leg S* = 1/2 lad- 
der that exists in the strong-coupling limit, J_l/J|| > 1. 
High-field, low-temperature magnetization, M{H < 30 
T, T > 0.7 K), data of single crystals and powder sam- 
ples have been fit to obtain J± — 13.3 K, J|| = 3.8 K, 
and A ~ 9.5 K, i.e. at the lowest temperatures finite 
magnetization appears at Hd = 6.6 T and saturation 
is achieved at Hc2 = 14.6 T. An unambiguous inflec- 
tion point in the magnetization, M{H,T = 0.7 K), and 
its derivative, dM/dH, is observed at half the satura- 
tion magnetization, Ms/2. This feature is symmetric 
about Ms/2, consistent with expectations for a simple 
spin ladder. Any presence of asymmetry, as was observed 
in Cu(Hp)Cl (^-|^, most likely arises from other factors. 
Our Ms/2 feature cannot be explained by the presence 
of additional exchange interactions, e.g. diagonal frus- 
tration Jp, but is well described by an effective XXZ 
chain, onto which the original spin ladder model (for 
strong-coupling) can be mapped in the gapless regime 
Hd < H < Hc2 [|2l|. After determining Hd and with 
no additional adjustable parameters, the magnetization 
data are observed to obey a universal scaling function . 
This observation supports our identification of BPCB as 
a two-leg S* = 1/2 Heisenberg spin ladder with J' ^ J||. 

The crystal structure of BPCB has been determined to 
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be monoclinic with stacked pairs oi S = 1/2 Cu^+ ions 
forming magnetic dimer units [^. The CuBr^^ tetra- 
hedra are co-crystalhzed along with the organic piperi- 
dinium cations so that the crystal structure resembles 
a two- leg ladder, Fig. 1. The rungs of the ladder are 
formed along the c*-axis (19.8° above the c-axis and 
+23.4° away from the a-c plane |^0|) by adjacent flat- 
tened CuBr^ ^ tetrahedra related by a center of inversion. 
The ladder extends along the a-axis with 6.934 A between 
spins on the same rung and 8.597 A between rungs. 
The magnetic exchange, J±, between the Cu^^ spins on 
the same rung is mediated by the orbital overlap of Br 
ions on adjacent Cu sites. The exchange along the legs 
of the ladder, Jy , is also mediated by somewhat longer 
non-bonding (Br- ■ -Br) contacts and possibly augmented 
by hydrogen bonds to the organic cations. A frustrating 
diagonal exchange, Jp, is possible, although it should be 
weak {Jp <^ J||), and so the potential of a finite Jp on 
the short diagonal was considered in our analysis. 




FIG. 1. Schematic of the crystal structure of BPCB. The 
legs (rungs) are along the a-axis (c*-axis), see text. 

Shiny, black crystals of BPCB were prepared by 
slow evaporation of solvent from a methanol solution of 
[(pipdH)Br] and CuBr2, and milling of the smallest crys- 
tals was used to produce the powder samples. The sto- 
chiometry was verified using CHN analysis, and 9 GHz 
ESR measurements were completely consistent with the 
previously reported data, i.e. ^(powder) = 2.13 |2^. In 
addition, deuterated specimens were produced and used 
in neutron scattering studies performed at the HFIR at 
Oak Ridge National Laboratory. No evidence for long 
range magnetic order or structural transitions was ob- 
served down to 11 K by powder diffraction and 1.5 K for 
single crystal diffraction in the [h 1] scattering plane. 
The low field magnetic measurements were performed us- 
ing a SQUID magnetometer. The high field work was con- 
ducted at the NHMFL using a 30 T, 33 mm bore resistive 
magnet and a vibrating (82 Hz) sample magnetometer 
equipped with a Cernox thermometer p3]. 



The low field, 0.1 T, magnetic susceptibility, x, of a 
powder sample, 166.7 mg, is shown as a function of tem- 
perature in Fig. 2. The data from single crystals, with 
the magnetic field oriented along the a-, b-, and c-axes in 
separate measurements, are indistinguishable from the 
results obtained with the powder specimen. The gen- 
eral shape of the curve is typical of low dimensional 
magnetic systems, and more specifically, it possesses a 
rounded peak at « 8 K and an exponential temperature 
dependence below the peak. Consistent with the neu- 
tron scattering results, no evidence of long-range order- 
ing was observed down to 2 K. A small extrinsic Curie- 
like impurity contribution (= 1.5% of the total number 
of Cu spins) and a temperature- independent diamagnetic 
term {xdia = —2.84 x 10~^ emu/mol, which is the sum 
of the core diamagnetism, estimated from Pascal's con- 
stants to be —2.64 x 10^'' emu/mol, and the background 
contribution of the sample holder) were subtracted from 
the data in Fig. 2. The Curie-Weiss temperature 6, 
and the Curie constant, C, can be extracted from a fit 
[x(T) = Xd^a + C/iT + 0), 50 K < T < 300 K], and we 
find C = 0.433 ± 0.002 emu K/mol and 6* = 5.3 ± 0.1 
K ||2^. These values are close to C = 0.425 emu K/mol 
(5 = 1/2, g = 2.13) and 6* « ( Jj_ + 2Jji)/4 ^ 5.2 K ||. 
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FIG. 2. The x{T) of a powder sample (166.7 mg) in 0.1 T. 
The line is the result of an exact diagonalization of a ladder 
Hamiltonian with 12 spins when J± — 13.3 K and Jy = 3.8 K, 
see text. The inset shows the M{H, T = 1 K) expectations of 
an exact diagonahzation of the alternating chain and ladder 
Hamiltonians with the exchange values given in the text. 

Initially, using exact diagonalization methods with 
12 spins, the x{T) data were fit to obtain the values 
Jj^ = 13.3 K, J|| = 3.8 K for a ladder Hamiltonian and 
J I = 13.7 K , J2 = 5.3 K for an alternating chain Hamil- 
tonian. Both fits are indistinguishable from the solid line 
shown in Fig. 2. Therefore, using only the low field 
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x{T) data, we were unable to distinguish between the 
ladder and alternating chain models, and this situation 
was not improved by fitting the M(i? < 5 T, T = 2 
K) data. However, in extensions up to the saturation 
magnetization, Mg, the alternating chain model gener- 
ated M{H,T < J||) curves that were asymmetric about 
as reported for Cu(Hp)Cl [||J2|], and the spin lad- 
der description predicted symmetric behavior, see Fig. 2 
inset. Since our experimental resolution was estimated to 
be sufficient to allow us to differentiate between the two 
models, the high magnetic field studies were initiated. 

The high field, < 30 T, magnetization of a powder 
sample, 208.2 mg, is shown in Fig. 3. Since Ms was 
reached in our studies, we were able to directly measure 
and subtract a small, temperature- independent contribu- 
tion {xdia ~ —2.84 X 10"'* emu/mole), which is the same 
value obtained in our low field work. Measurements were 
also made on a single crystal, 18.9 mg, with H \\ a-axis 
and for T > 1.6 K. Within the resolution, the data are 
the same for the powder and single crystal samples. Fur- 
thermore, the data were acquired while ramping the field 
in both directions, and no hysteresis was observed. 

The low energy states of the spin ladder Hamiltonian 
can be mapped, in the strong-coupling limit, onto the 
5 = 1/2 XXZ chain H], allowing M{H,T) to be mod- 
eled. The solid lines in Fig. 3 were obtained by numerical 
integration of the Bethe ansatz equations for the effec- 
tive XXZ chain , using the parameters describing the 
spin ladder fit for x(T). All of the data, Figs. 2 and 3, 
are reproduced by one set of exchange values when using 
the ladder model. On the other hand, the alternating 
chain model fails to fit all of the data with a single set 
of parameters. For example, the dashed line in Fig. 3 
is M/Ms{H,T = 0.7 K) calculated from the alternating 
chain mapping onto the XXZ chain when Ji = 13.7 K 
and J2 = 5.3 K, i.e. the values obtain from fitting x(^)i 
Fig. 2, by an alternating chain model. In addition, our 
data were analyzed with a ladder model that also incor- 
porated a frustrating interaction, Jp p6| , and we can 
estimate an upper bound of Jp < 0.5 K. Consequently, 
all of the data are consistent with a strongly coupled lad- 
der description for BPCB, where J± — 13.3 ± 0.2 K, and 
J|| = 3.8 ± 0.1 K. 

To leading order, q^bHci ^ J± — J\\, and gfiBHc2 — 
J± + 2 J|| [p6| , p7| . Using the previously mentioned param- 
eters, we obtain Hd = 6.6 T and Hc2 = 14.6 T, identical 
with the experimental results. The inset in Fig. 3 shows 
the derivative curve, d{M/Ms)/dH, of our data at the 
lowest temperature. The symmetric double bump struc- 
ture and its evolution with temperature has been studied 
theoretically [T^ but has not been observed previously in 
S" = 1/2 two- leg ladder materials. Even though our the- 
oretical curve somewhat overestimates the sharpness of 
d(M/Ms) /dH, the overall agreement between theory and 
experiment, including the evolution of M(H,T), Fig. 3, 
is excellent, and involves no adjustable parameters once 



Hci is defined. Furthermore, the fact that we see only 
one feature at Ms/2 between Hd and Hc2 is evidence 
that our strongly interacting dimers are not coupling to 
form 2D ||2il or 3D ||29|,|30|l networks. 
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FIG. 3. The normalized magnetization, M/M^, of a pow- 
der sample (208.2 mg). The data traces are limited to « 150 
(of ~ 3000) points for clarity. The lines are spin ladder pre- 
dictions of an effective XXZ chain when J± — 13.3 K and 
J|| = 3.8 K. At T = 0.7 K, the inflection point at Ms/2 is 
clearly visible, and the inset shows the derivative of this data. 
The dashed line is the alternating chain model prediction for 
0.7 K when Ji, J2 are taken to be the values obtained from 
fitting the data in Fig. 2, see text. 

At Hd, BPCB undergoes a transition from gapped 
dimer pairs to a gapless Luttinger liquid phase with 
fcrmionic excitations, where the magnetization is pro- 
portional to the fermion density P, ^ , ^ . This transi- 
tion can be described as a condensation of a dilute gas of 
bosons (dimers) , and quasiparticle interactions are irrele- 
vant at the transition point. At Hc2, an analogous situa- 
tion exists where the transition is between the Luttinger 
liquid and spin polarized phases. When T, g^slH — Hd\, 
and gnB\Hc2 ~ H\ are ^ Jy, the ID magnetization is 
predicted to obey the universal scaling law (assuming 
J±/J\\ ^ 1) that may be written as: 



M{H,T) 
Ms 



2fcBr/J|| M{g^iB[H - Hd]/kBT) , 



M{H,T) 
Ms 



2kBT/.h M{g^iB[Hc2 - H]/kBT) 



where the universal function M. is the fermion density . 
This theoretically predicted scaling behavior is compared 
to the data in Fig. 4, where the agreement is impressive. 
It is important to stress that the scaling shown in Fig. 
4 has been theoretically predicted |^ and is not a result 
of extracting scaling variables on the basis of the data 
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p2[ . In an isolated spin ladder, scaling is expected at the 
lowest temperatures, T ^ J\\. A deviation from scaling is 
observed for T = 0.7 K, which suggests that other weak 
interactions, such as Jp or J', may begin to have a subtle 
influence, while the data up to 4.47 K appear to obey the 
scaling theory. The T^/"^ scaling of the magnetization at 
the critical point H = Hd, Fig. 4, is further evidence 
that BPCB is a two- leg spin ladder with J' < J|| Q. 
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FIG. 4. The scaled data in the vicinity of Hci and Hc2- 
The solid lines are the predictions of the theory when Hci is 
fixed. The inset shows the T^^^ scaling behavior at Hd- 

In summary, analysis of M{H < 30 T, T > 0.7 K) has 
allowed us to identify BPCB as a two- leg 5=1/2 spin 
ladder in the strong-coupling limit, Jj_/J\\ ~ 3.5. A sin- 
gle set of exchange constants, J± = 13.3 K and J|| — 3.8 
K, are able to accurately describe all of the data. The 
M{H « Hci or Hc2,l K < T < 4.5 K) data exhibit 
scaling behavior in the universality class of the ID dilute 
Bose gas transition Although we have considered 

the potential existence of additional exchange interac- 
tions Jf and J', effects arising from these parameters 
are not prominent in the present data. However, since 
subtle differences arise between the theoretical predic- 
tions and the data at the lowest temperature, additional 
perturbing interactions may be present. 
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